A supercritical series analysis for the generalized 
contact process with diffusion 

W.G. Dantas 

Institute de Fisica, Universidade de Sao Paulo, Caixa Postal 66318 05315-970 Sao 
Paulo, Sao Paulo, Brazil 

J.F. Stilck 

Instituto dc Fisica, Universidade Federal Fluminense, 24210-340, Niteroi, RJ, Brazil 

Abstract. We study a model that generalizes the CP with diffusion. An additional 
transition is included in the model so that at a particular point of its phase diagram 
a crossover from the directed percolation to the compact directed percolation class 
will happen. We are particularly interested in the effect of diffusion on the properties 
of the crossover between the universality classes. To address this point, we develop 
a supercritical series expansion for the ultimate survival probability and analyse 
this series using d-log Pade and partial differential approximants. We also obtain 
approximate solutions in the one- and two-site dynamical mean-field approximations. 
We find evidences that, at variance to what happens in mean-field approximations, 
the crossover exponent remains close to </) = 2 even for quite high diffusion rates, and 
therefore the critical line in the neighborhood of the multicritical point apparently does 
not reproduce the mean-field result (which leads to (j) — 0) as the diffusion rate grows 
without bound. 
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1. Introduction 

Phase transitions in stochastic models have attracted a great attention in recent 
years. Ahhough much work has been done on systems whose stationary states are 
described by thermodynamics, also phase transitions in systems far from equilibrium 
have been studied in much detail, and it is remarkable that concepts of scaling and 
universality, developed in the context of thermodynamic phase transitions, have been 
applied successfully in these new situations as well. In particular, systems which exhibit 
absorbing states, and therefore do not obey detailed balance, are a rich group of models 
where out of equilibrium phase transitions have been studied. Much of this work has 
been done based on numerical simulations, stochastic models are very well fitted to 
these techniques, but analytical approaches are also useful and may sometimes lead to 
very precise results. Therefore, as in other fields of Physics, simulational and analytical 
approaches are complementary to study non-equilibrium phase transitions [1]. 

One of the largest and most studied universality classes in models with absorbing 
states is the directed percolation (DP) class, and in this class the most studied model is 
the contact process (CP), which in one of his interpretations was formulated as a very 
simple model to describe the evolution of an epidemic disease, which spreads through 
contact between healthy and sick individuals, placed on sites of a lattice. In this model, 
at each site of the lattice an individual is located, and this individual will be in one 
of two states: healthy or sick. It is usual to associate a sick individual to a particle 
and a healthy one to a hole. Particles are created autocatalytically, that is, a particle 
may be created at an empty site with a rate which is proportional to the number of 
occupied first neighbor sites of this site. An occupied site may become empty with a 
unitary rate (spontaneous healing). Even in one dimension, a phase transition is found 
in the model between the absorbing stationary state (in which all individuals are healthy 
and no new sick person may be generated by contact, or else for a lattice without any 
particle) and a stationary phase with a nonzero fraction of particles, called active state 

The contact process is related to several other stochastic models [HI HJ |5l |6] , and 
the directed percolation universality class seems to include all models where a transition 
between an absorbing and an active state occurs, with a scalar order parameter, short 
range interactions and no conservation laws [7j. No exceptions to this conjecture have 
been reported so far [8]. 

Recently, a generalization of the CP was studied using mean-field approximations 
and simulations [9], as well as series expansions [TO]. In this model, an additional 
process is included besides the autocatalytic creation and the spontaneous annihilation 
of particles: the autocatalytic creation of holes, that is, an occupied site may become 
empty by two processes: either spontaneously or with a rate which is proportional to 
the number of empty first neighbor sites. It is more natural to describe this model in 
terms of two distinct types of particles, A and B, requiring that each lattice site at any 
time is either occupied by an A or a i? particle. The process B ^ A occurs with a rate 
which is proportional to the number of first neighbor sites with A particles, whereas the 
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opposite process A ^ B may happen either spontaneously or with a rate proportional 
to the number of first neighbor sites occupied by B particles. The CP is a particular 
case of this model where no autocatalytic creation of B particles is allowed. In another 
particular case, an additional symmetry is present in the model: when the spontaneous 
creation of B particles is suppressed, the model becomes symmetric with respect to 
the interchange of A and B particles for equal creation rates of A and B particles, 
being known as the biased voter model [11], with a transition between two symmetric 
absorbing states, where the lattice is totally filled with the same type of particle. If the 
rates of the two processes {A B and B ^ A) are the same, the density of A and B 
particles does not change as the system evolves, but if the rate of one of these processes 
is larger, the system will reach the absorbing state in which only the particles created 
at a larger rate are present. This model, having an additional symmetry, belongs to the 
compact directed percolation (CDP) class, with critical exponents which are different 
from the ones found in the CP. The crossover between both universality classes was 
the major motivation to study this generalized model. In particular, the results of a 
two-variable series expansion for the ultimate survival probability of the model lead to 
rather precise evidences that the crossover exponent at the point where the universality 
class of the model changes is given by = 2, coincident with the mean-field value, 
which agrees both with the limits obtained by Liggett |12] and with simulational results 
[T3] for the same crossover in the Domany-Kinzel automaton [H], thus providing an 
evidence that the process of update (parallel in the DK automaton and sequential in 
the CP) does not change the value of the crossover exponent. 

It is of interest to study the effect of diffusion on the behavior of such models. 
This was already done for the CP, with the result that the transition point becomes a 
critical line as the new variable (rate of diffusion) is introduced in the model IT6] . In 
this C3iSG, clS expected, mean-field results are obtained in the limit where the evolution 
of the system is dominated by diffusion, both in the values of the transition rate and 
of critical exponents. Thus the limit of infinite diffusion rate may also be viewed as 
a crossover transition between non-classical and classical behavior. These results may 
be heuristically justified if we note that in this limit, since diffusion processes are much 
more probable that the others, for each other process the local densities may be replaced 
by the global ones, eliminating the effect of fluctuations. 

In this work we include diffusion in the generalized CP model described above, 
in order to find out its effect on the phase diagram and critical exponents of the 
model. After defining the model more precisely in section [21 we obtain its phase 
diagram in the two-site cluster mean-field approximation in section [3l We then proceed 
obtaining supercritical series expansions for the survival probability in section [U which 
are analysed using Pade approximants and two- variable partial differential approximants 
in section [5l Conclusions and final comments may be found in section [6] 
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2. Definition of the model 



The model is defined on a one-dimensional lattice with N sites and periodic boundary 
conditions. Each site is occupied either by a particle A or a particle B, no holes are 
allowed. The microscopic state of the model may thus be described by the set of binary 
variables rj — (771, 772, ■ ■ ■ , Vn), where 771 = or 1 if site i is occupied by particles B or A, 
respectively. 

The model evolves in time according to the following Markovian rules: 

(i) A site i of the lattice is chosen at random. 

(ii) If the site is occupied by a particle B, it becomes occupied by a particle A with a 
transition rate equal to PaTiA/'^, where tia is the number of A particles in the sites 
which are first neighbors to site i. 

Pa/2 



ABB 



A A B 



(iii) If site i is occupied by a particle A, it may become occupied by a particle B through 
three processes: 

• Spontaneously, with a transition rate Pc- 



A 



Pc 



B 



Through an autocatalytic reaction, with a ratep^^e/S, where 7^B is the number 
of B particles in the sites which are first neighbors to site i. 

Pb/2 



A A B 



ABB 



By interchanging the B particle with an A particle located in a first neighbor 
site to site i. This transition occurs with rate equal D. 



A 



B 



D 



B A 



where □ indicates the state of the first neighbour. 

We define the time in such a way the the non-negative rates Pa, Pb, and Pc obey the 
normalization Pa + Pb + Pc ~ 1- We may then discuss the behavior of the model in the 
{Pa-iPc-i D) space without loss of generality. 

When Pc 7^ the configuration where all the sites are occupied by particles B 
is absorbing, while for p^ — both symmetrical configurations in which all sites are 
occupied by the same particles are absorbing, independently of the value of D. In the 
particular case pf, = 0, in which the model corresponds to the CP with diffusion, there 
exists a critical line in the {pa,Pb, D) space where an active configuration (with nonzero 
density of A particles) becomes identical with the absorbing one. This line becomes a 
critical surface in the general case p^ 7^ 0. Based on the symmetries of the model this 
surface must be of second order with critical exponents of DP class for finite values 
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of the diffusion rate D, except for the case Pc = 0. In this case, the transition occurs 
between two absorbing states (lattice full by particles A or B) at the the point Pa = 1/2, 
\fD, with critical exponents of the CDP class. 

For a certain fixed value of D, the behavior of any stationary density close to 
multicritical point {pa = 1/2, pc = 0) should exhibit the scaling form. 



9(p. - 1/2, p„ D) ~ (p., - iftY-'fi^F ■ (1) 

The critical exponent associated with the density variable g, eg{D), should 
correspond to the CDP universality class, and the scaling function F{z) is singular 
at a value Zq{D) of its argument, which corresponds to the critical line for a given value 
of D. Thus, the critical line is asymptotically given by pc = zo{D){pa — 1/2Y^^\ In 
this way we have two exponents as a function of the diffusion and it is interesting to 
find out how the exponents change as diffusion processes are introduced in the model. 



3. Cluster Approximations 

We have derived solutions for the cluster dynamic approximations in simple mean-field 
(one-site) and two-site levels [1]. The simple mean-field solution is independent of the 
diffusion rate, and the critical line is located at = 1/2, independently of Pc- Thus, 
in this approximation the crossover exponent vanishes identically. The lowest order of 
mean-field cluster approximation in which the effect of diffusion is present in the results 
is the two-site approximation. Without going into the details of this approximation, 
since the calculations are similar to the ones performed recently in the model without 
diffusion [9], in the two-site level cluster approximation the critical line is given by: 



Pa 



4 - (pe + 2D) + V(Pc + 2Df + 8p, 



(2) 



We notice that as the value of D grows, the critical line approaches the result 
obtained in the one-site approximation, as expected. If we expand equation [2] for small 
values of Pc, we obtain — 1/2 f» 2pc/ D, showing that the two-site approximation leads 
to the crossover exponent = 1 for any finite value of the diffusion rate D and = in 
the limit of infinite diffusion rate. For vanishing diffusion rate, we obtain the crossover 
exponent = 2. 

Another point which is worth noticing is the crossover as the diffusion rate 
approaches infinity, where mean-field behavior is expected. For this purpose, we may 
define the variables a = {1—Pa)/Pa and -De// = aD/{l + aD). Rewriting the expression 
for the critical condition [2] in terms of the variables x=l — a,y = l — D^jf, and pb, we 
obtain: 

- 2) 

PbX^ + (1 - 4:pb)x^ + {5pb -l)x+l-2pb ^ ^ 
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It is apparent that, as long as pi, < 1/2, we have 0=1 for the crossover between the 
regimes of finite and infinite diffusion. In the particular case pf, = 0, which corresponds 
to the CP, we recover two-site approximation results obtained before [Mj. For the 
particular case pf, = 1/2, which corresponds to the CDP limit, we have the solution 
X = for any value of y, so that the locus of the transition is not affected by the 
diffusion. These results are illustrated in figure [H 




Figure 1. Results of the two-site approximation. In the left panel we see different 
phase diagrams for fixed values of pb using the new variables a and Def / . In the right 
we have the estimates of the crossover exponent in each of these diagrams. In all cases 
the parallel dashed lines indicate that this exponent is (j) — I for any value of the 
parameter ph < 1/2. 



4. Derivation of the supercritical series for the model 

Now let us develop a three-variable supercritical series expansion for the model. We 
follow closely the operator formalism presented in the paper by Jensen and Dickman on 
series for the CP process and related models [I7j. We may represent the microscopic 
configurations of the lattice by the direct product of kets 

\v) = <^h), (4) 

i 

which are defined to be orthonormal 

{vW) = ]lK,v.'- (5) 
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Now we may define A particles creation and annihilation operators for the site i: 
Al\rji)^{l-rii)\r]i + l), 

AilVi) = VilVi - 1)- (6) 
In this formalism, the state of the system at time t may be represented as 

\m) = T.piv,t)\v)- (7) 

If we define the projection onto all possible states as 

(I^E(^h (8) 

{ri} 

the normalization of the state of the system may be expressed as { l^p) = 1- In this 
notation, the master equation for the evolution of the state of the system is 

'^ = sim- (9) 

The evolution operator S may be expressed in terms of the creation and annihilation 
operators as 5" = XR + V where 

R = E[7(2 - 4-1 A-i + + i]{Ai - 4a) + 

i 

+ D ^(1 - 4_i A)^^-l4 + (1 - 4+1^0 A+i4, (10) 

i 

V = E(4 + 4a - i)(4-iA-i + 4+i^.+i), (11) 

i 

and the new parameters 

Pa 
Pb 

2pc 



7 



and 



A 

were introduced. 

We notice that the operator R includes the diffusion and the annihilation of A 
particles processes, while the creation of A particles is present in the operator V. 
Therefore, at small values of the parameter A creation of A particles is favored, and 
the decomposition above is convenient for a supercritical perturbation expansion. Let 
us show explicitly the effect of each operator on a generic configuration (C) . 

m = 7 E(c:) + 27 Hic'j) + E(cD + 



+ D 



E(C) + E(cr) + E(C') 



+ 



- b{ri + 2r2) + D{t, + 2^2) + r](C), (12) 
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where the first sum is over the vi sites with A particles and one B neighbor, the second 
sum is over the r2 sites with A particles and two B neighbors, the third sum is over all 
r sites with A particles of the configuration (C). The next sum is again over the ri sites 
with A particles and one B neighbor and the two last sums are over the r2 sites with A 
particles and two B neighbors. Configuration (C-) is obtained replacing the A particle 
at site 2 by a B particle, (Cf ) is obtained interchanging the A particle at the site i with 
its single B neighbor and finally, {C"'^"^'^^) is a configuration where the A particle at the 
site i is interchanged with the B particle located at the right (r) or to the left (/) of site 
i. It is convenient associate the diffusion with the annihilation process to avoid some 
ambiguity in truncating the series in a certain order. In the other hand, the action of 
operator V is 



(13) 



where the first sum is over the qi sites with B particles and one A neighbor, the second 
sum is over the 52 sites with B particles and two A neighbors. Configuration (Cf ) is 
obtained replacing the B particle at site i in configuration (C) by a A particle. 

Since we are interested in the long time behavior of the system, it is useful to 
consider Laplace transforms. For the state ket we have: 



(14) 



and inserting the formal solution \'ip{t)) = e \iIj{0)) of the master equation [9l we find 

ms)) = {s-s)-'\m)- (15) 

The stationary state lip^oo)) = lim^^^oo |^(i)) may then be found noticing that 

1^^(00)) = lims|^(s)), (16) 

which may be obtained integrating [H] by parts. A perturbative expansion may be 
obtained assuming that \ip{s)) may be expanded in powers of A and using | 



Since 



we arrive at 



|^(s)) = |^o) + A|^i) + A2|^2) + 



1 + A 



V -\R 



\m)- 



[17) 



s-V -XR s-V 



s-V 



R + X' 



:R' + -- 



1^1) 
1^2) 



s-V 
1 

s-V 
1 

s-V 



\m) 

R\^o) 



(19) 



(20) 
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The action of the operator {s — V) ^ on an arbitrary configuration (C) may be found 
noting that 

(. - V)-\C) = s-\C) + ^^^^(C), (21) 
and using the expression [T3l for the action of the operator V, we get 



v)-\c) = sJic) + {s-vy 



E(cr) + 2E(c;) 



(22) 



where the first sum is over the qi sites with B particles and one A neighbor, the second 
sum is over the q2 sites with B particles and two A neighbors, and we define Sg = 1/ (s+g), 
where q = qi + 2g2- 

It is convenient to adopt as the initial configuration a translational invariant one 
with a single A particle (periodic boundary conditions are chosen). Now we may notice 
in the recursive expression [22] that the operator {s — V)^'^ acting on any configuration 
generates an infinite set of configurations, and thus we are unable to calculate {ip) 
in a closed form. We may, however, calculate the extinction probability p{s), which 
corresponds to the coefficient of the vacuum state |0). As happens also for the 
models related to the CP studied in [T7| configurations with more than j particles 
only contribute at orders higher than j, and since we are interested in the ultimate 
survival probability for A particles Poo = 1 — lims^osp(s), Sg may be replaced by 1/g 
in Eq. [22l An illustration of this procedure may be seen in [TO] . 

The algebraic operations above may be easily performed in a computer using 
a proper algorithm. The configurations are expressed as binary numbers and the 
coefficients as double precision variables. With rather modest computational resources 
(Athlon MP2200, double processor, 1Gb memory) it is not difficult to calculate the 
coefficients up to order 22. The required processing time amounts to about 2 hours, 
the limiting factor is actually the memory required for the calculation. We define the 
coefficients bij^k as: 

-1-1 22 i-l i-2 

Poo = 1 - -(27 + 1)A - -(27 + - E E E KA'i'DK (23) 

^ ^ i=3 j=0 k=0 

The case a = is leads to a series which is identical the the one for the CP with diffusion 
obtained in [15]. Also, for Z) = the two- variable series for the model without diffusion 
is |10j is recovered. The set of coefficients is too large to be included here, but may be 
obtained upon request from the authors. 



5. Analysis of the series 

To obtain estimates of the critical properties of this model from the supercritical series 
for the ultimate survival probability, given by equation [231 we start by using the d- 
log Pad approximants approach. These approximants are defined as ratios of two 
polynomials 
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T 
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iVi 


Ac 


U 


U.z 
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U.4y4 1 / 






ii 
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1 
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^047fi 






11 


11 


0.50476 


0.7 


0.2 


10 


10 


0.64708 






11 


11 


0.64750 


0.8 


0.2 


10 


10 


0.75146 






11 


11 


0.71941 



Table 1. Examples of estimates for critical points obtained by d-log Pad approximants. 
Notice that as the value oi D — D/{1 + D) grows the dispersion in the estimates 
obtained from different approximants for this critical point also grows. 



FlmW = = . = fW- (24) 



In our case the function /(A) represents the series for ^ In Poo (A), where we fix the two 
remaining parameters 7 and D. Therefore, we may obtain approximants with L + M < 
22, and since it is known that diagonal (L = M) and near-diagonal approximants 
usually exhibit better convergence properties, we restrict our calculations to the set of 
approximants such that L = M+6, with 6 = 0, ±1. We built approximants for the values 
of 7 ranging between and 45, using D = D/{1 + D) = 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.7, 0.8 
and 0.85. In the table [1] some estimates for the critical value of the parameter A are 
shown. The estimates poles are roots of the polynomial Qm (poles of the approximant). 

With the data provided by the set of approximants, we obtained the critical lines 
using the original variables pa and pc, for different values of the diffusion rate D. Each 
point of these curves, shown in the figure [21 was calculated as the average of the estimates 
provided by the set of approximants, and the error bar associated to it corresponds to 
the standard deviation of the estimates. As the diffusion rate increases, the dispersion of 
the approximants also grows and the precision of the result becomes smaller generating 
relatively higher errors. Even so, it is possible to notice that all the critical curves 
approach of the multicritical point with a quadratic curvature, and thus the present 
results are different from the ones provided in the mean-field two-sites approximation 
shown above, which lead to = 1. 

The critical exponent f3 associated to the order parameter also can be estimated 
from the d-log Pad approximants, calculating the residue associated to the physical 
pole of each approximant. In the figure [H the average values of the estimates for /3 
are depicted for different values of the diffusion rate. They are consistent with the 
conjecture that the value of the DP universality class Pdp = 0.276486(8) [H] applies to 
the model with finite rates of diffusion. Also, as expected, we notice that the dispersion 
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□-€]D = 0.1 
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^^D = 0.3 




^TcD = 0.5 
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D = 0.8 




'v-vD = 0.85 





0.3 



0.4 




Figure 2. It the top panel we see the phase diagrams of the model obtained from d-log 
Pad approximants for several values of the diffusion rate D = 5/(1 + Z)). The graph 
in the bottom shows log-log plot of the same curves with the dashed line indicating 
the value of (/)~^ = 1/2. 



of the estimates grows as the CDP point is approached (pc --^ 0) 
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Figure 3. Behavior of the exponent /? as a function of the parameter pc for three 
different values of the diffusion rate. The value corresponding to the DP universality 
class is shown as a dashed line. 



The behavior of the model in the limit of infinite diffusion rate is out of the reach 
of this analysis, partially because the reduction of the series to a variable leads to poor 
results in the neighborhood of a multicritical point, as is known of other similar cases 
[To] and therefore the d-log Pad approximants suffer great fluctuations in the region of 
high diffusion, as already found in the model without diffusion |T0]. 

In an attempt to improve the series analysis we may use the Partial Differential 
Approximants (PDA' s), which are appropriate for the study of two- variables series 
of models which undergo multicritical phenomena [19]. They may be regarded as a 
generalization to two variables of the d-log Pade approximants. The defining equation 
of the approximants is 

P^ix, y)Fix, y) = Q^{x, v)^^^^ + ^n(x, 1/)^^^, (25) 

where P, Q, and R are polynomials in the variables x and y with the set of nonzero 
coefficients L, M, and N, respectively. The coefficients of the polynomials are obtained 
through substitution of the series expansion for the quantity which is going to be 
analyzed 

f{x,y)= fik,k')x'y'' (26) 

k,k'=0 

into the defining equation [25] and requiring the equality to hold for a set of indexes 
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defined as K. The coefficients of tlie polynomials are then found solving a set of linear 
equations. Since the coefficients fk^k' of the series are known for a finite set of indexes 
this sets an upper limit to the number of coefficients in the polynomials. The number of 
equations has to match the number of unknown coefficients, and thus K = L + M+N — 1 
(as in the Pade approximants, one coefficient is fixed arbitrarily). An additional issue, 
which is not present in the one-variable case, is the symmetry of the polynomials. Two 
frequently used options are the triangular and the rectangular arrays of coefficients. The 
choice of these symmetries is related to the symmetry of the series itself |20j . 

Let us suppose that the quantity represented by the series is expected to have a 
multicritical behavior at a point {xcVc), described by the scaling behavior: 

f{x,y)^\ASrfZ^l^y (27) 

where 

Ax = {x- Xc) - {y- i/c)/e2, 

and 

Ay = {y- Vc) - ei (x - Xc) . (29) 

Here e/ is the critical exponent of the quantity described by / when = 0, ei and 62 
are the scaling slopes [19] and is the crossover exponent. The function Z{z) is singular 
for one or more values of its argument, corresponding to the critical line(s) incident on 
the multicritical point. Once the coefficients of the defining polynomials are obtained, 
the estimated location of the multicritical point corresponds to the common zero of the 
polynomials Qm and Rjy. This may be seen substituting the scaling form [271 in the 
defining equation [25] of the approximant. The exponents and scaling slopes may also 
be obtained directly from the polynomials, without integrating the partial differential 
equation. A detailed discussion of the algorithm, as well as computer codes, may be 
found in [20]- 

Since in our original series we have three variables, to accomplish the analysis using 
PDA's, we fix the value of the variable D and we generate series in the following the 
variables x = X and y = jX, to avoid numerical errors due to the divergence of the 
variable 7 in the multicritical point of crossover to the GDP universality class. With 
that, for each value of D we calculate about 22 approximants, choosing different sizes 
and configurations for the set {K, L, M, N}, but maintaining the constraint K ?s M. 
The results, calculated as averages over the estimates of the set of approximants, for 
the exponents Cf = (3' and as functions of the diffusion rate D are shown in the 
figure [H The same quantities were calculated imposing as a constraints the location 
of the multicritical point, fixed at x = and y = 1. The values of the exponents 
with this constraint are represented as squares of the figure [H and are very close to the 
non-constrained estimates (represented by circles). Again the error bars are estimated 
from the dispersion of the results in the set of approximants. The crossover exponent (f) 
seems to be invariant with the diffusion rate, as already be suggested by the results of 



(28) 
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the d-log Pad approximants. Actually, the greatest deviations from = 2 are seen for 
intermediate values D = D/{l+D) ~ 0.5. In fact, for D = 0.1 we have = 2.002±0.051 
{(f) = 2.038 ± 0.201, for the constrained case) while for D = 0.5, = 2.134 ± 0.415 
(0 = 2.105 ± 0.201, for the constrained case) and in D = 1.0 we have = 2.025 ± 0.090 
(0 = 2.031 ±0.121). Surprisingly, in the region close to infinite diffusion rate, ot D = 1, 
the estimates of the crossover exponent show a dispersion of the same order of the one 
found in the region of vanishing diffusion D ^ 0. These result suggests that the crossover 
does not coincide with the mean-field result for any non-zero value of the diffusion rate. 
We notice that the value = 2 is inside the error bars for all the estimates. 
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Figure 4. Behavior of the exponents cj) and (3' as a function of the diffusion rate. The 
circles represents the non-constrained results and the squares denote the constrained 
results. The dotted line is only a guide for the eyes, while the dashed lines correspond 
to the values (p ^ 2 and /3' = 1 of the case without diffusion. 



Using the method of characteristics, the equation [251 defining the approximants may 
be integrated. A time-like variable r is introduced, and the partial differential equation 
will be equivalent to a set of two ordinary differential equations: 

^ = QM(a;(r),?/(r)), 
ar 

^ = R^(xiT),y{T)). (30) 

The characteristics are the trajectories (x(r),?/(r)) obtained solving these equations, 
and the estimate of the approximant on each characteristic may be found through 
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integration, once we know the value on a particular point. It is possible to show that any 
critical line which is incident on the multicritical point will itself be a characteristic curve, 
and we will use this result to estimate the DP critical line in the model. The critical 
lines obtained by the method of characteristics are shown in the figure [5], calculated for 
D = 0, 0.1, 0.2, 0.3, 0.5, 0.7, and D = 0.8. 
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Figure 5. Critical lines obtained by the method of characteristics for different values 
of the diffusion rate D. This approach allows us to estimate the curves close to the 
multicritical point of DP-CDP crossover. 



Once again, the approach fails to estimate the behavior for the limit of infinite 
diffusion rate. In fact, only values to up to D = 0.8 could be studied with reasonable 
precision. On the other hand, in [in], it is shown that a good estimate of the critical 
line is obtained by the scaling function x = zq{1 — y)'^. A similar procedure may be 
applied to the model with diffusion, adopting the Ansatz x = zq{D){1 — y)'^^^\ where 
0(D) = 2, VD and zo{D) is a parameter chosen to reproduce the result in a region were 
we have good estimates. In the figure [6] we show estimates to the critical lines from the 
method of characteristics and extended scaling curves for the diffusion rates D = 0.2, 0.5 
and D = 0.8. In the same figure, the extended scaling curves are depicted up to the limit 
of infinite diffusion rate. In that case, the curve has a maximum, unlike of the behavior 
of the other critical lines. Therefore, if the hypothesis that all the critical curves can 
be approximated by a scaling function such as x = Zq{D)(1 — yY^^^i is true, then in the 
limit of infinite diffusion rate only the limiting cases P6 = (that equals to the CP) and 
Pc = (voter model) have the same behavior that is predicted by the mean-field for this 
regime. 
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Figure 6. Top: Estimates of critical lines obtained by the method of characteristics 
(dashed lines) and using the extended scaling form (symbols) are compared. Bottom: 
Critical lines calculated using the extended scaling form x{D) = zq{D){1 — y)^ for 
different values of diffusion, including the infinite diffusion rate limit. 



To test if the hypothesis that the critical curves can be approximated by the function 
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x{D) = zo{D){l —yY'y^) for all value of the diffusion rate is reasonable, we calculate the 
behavior of -De// as a function of a, in the same way as was done to obtain the curves 
shown in the figure [H Remembering that a = (1 —pa)/Pa and -De// = a-D/(l + aD) we 
fix the value of the parameter ph and using the values of Zo{D) our results estimate the 
location the multicritical point {ac,Deff = 1). If mean- field behavior would be correct 
for all values of pb in the infinite rate of diffusion limit, then Oc should always be equal 
to one. However, as may be seen in figure U\ this value changes. We also notice that in 
the figure [7] an agreement between the points obtained of the critical curves estimated 
by the approach of the characteristics and the curves approached, giving support to the 
possibility of that ac varies with pf, and, consequently with the parameter pc, as shown 
in the figure [61 
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Figure 7. Phase diagrams using a and D^/f as variables. For each value of parameter 
Pb we have a different value for ac, that is the point where the curve achieves the infinite 
diffusion limit. Again the symbols are results from the method of characteristics and 
dashed lines are extended scaling functions. 



The curves for -De// as a function of a provide another evidence of the peculiar 
behavior of the critical lines in the limit of infinite diffusion rate. For the CP with 
diffusion, at the DP-mean field change of universality class, the estimated crossover 
exponent is = 4, as found in [16]. In the figure [8] we show as this value of the 
exponent is obtained using the same values for ac found in the right panel of the 
figure [3 

In summary, the results obtained by the PDA's calculated for distinct values of D 
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Figure 8. Log- log plots for critical curves with different values of pb, showing that 
the choices of ac from the results exhibited in figure [7| result in a crossover exponent 
^ = 4, as expected for the CP. 



and presented in the figure H] are consistent with what is known in the hterature [16] of 
the behavior expected for the crossover between the DP universahty and the mean-field 
behavior. The results presented here seem to support the conjecture that in the limit 
of infinite diffusion the critical lines of the model are coincident with the simple mean 
field result only in the extrema, which correspond to the usual CP and the voter model. 



6. Conclusion 

In the generalized model without diffusion, studied using series expansions in [10], it was 
found that the DP-CDP crossover exponent was very close to the two-site mean-field 
approximation result = 2. Here we found that when diffusion is introduced in the 
model, the same approximation leads to = 1 for any nonzero finite diffusion rate and 
= in the limit of infinite diffusion rate. As usual, this latter result is coincident with 
the one which is found applying a simple one-site mean-field approximation to the model. 
However, the series analyses for the model presented here support the conclusion that 
the introduction of diffusion does not change the crossover exponent. If this conclusion 
remains true in the limit of infinite diffusion, as some results above suggest, the critical 
line joining the points corresponding to the CP model and the voter model could not be 
the curve pa = 1/2 predicted by the one-site mean field approximation, since it should 
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be a quadratic curve in the neighborhood of the multicritical point at Pc = 0. 

It is well known that the Pade approximants will provide poor estimates of the 
critical parameters in the neighborhood of a multicritical point. Actually, this was one of 
the motivation for the development of PDA's, which are suited to estimate multicritical 
behavior [T9j. As expected, Pade approximants show an increasing dispersion of the 
estimates as the DP-CDP multicritical point is approached, whereas PDA's lead to 
good estimates in this region [10]. As diffusion is introduced in the model, we found out 
that PDA's apparently are reliable for low diffusion rates, but as the rates are increased 
again the dispersion of estimates provided by different approximants grows and for 
infinite diffusion rate and no longer obtain reliable estimates from the approximants. 
One may suppose that, similar to the poor performance of Pade approximants close 
to multicritical points, the PDA's also fail as the multicritical point of higher order is 
approached for infinite diffusion rate. A generalization of the PDA's to handle a three 
variable series may be helpful to study this limit, and we are presently working in this 
direction. 

Finally, in our opinion the conclusion that the crossover exponent does not change, 
even in the limit of infinite diffusion rate, which our results seem to support, should be 
viewed with some caution. We stress that in order to obtain a definite series expansion 
for the model, as was also necessary in earlier studies of series expansions for similar 
models with diffusion, we were forced to include the diffusive term into the part of 
the evolution operator which is treated as a perturbation, and therefore it might be 
possible that conclusions for high diffusion rates are misleading. Nevertheless, for low 
diffusion rates we may have more confidence in the results from the analysis of the series 
expansion, and there are clear evidences that, unlike what happens in the two-site mean 
field approximation, the value 0^2 found in the absence of diffusion, is still valid when 
diffusive processes are allowed. 
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